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‘ ABSTRACT 

The objotive of this study was to select the optimal combination of 
threshold values for tho 3x3 pixel spatial filtering for post-classification 
implnntod in Image- 100 system of INPE and to evaluate the effects of this 
procedure on the accuracy of wheat area wliich was estimated by Image-lOO 
system using a hybrid classifier and LANDSAT digital data obtained from a 
single pass. An area of 800 hnf in Cruz Alto which is ono of the most important 
municipals for wheat production in southern Brazil,, was selected for this 
study. Different threshold combinations; (1,2), (2,10, (2,2), (2,3), (3,2), 
(3,4) and (3,5) were employed in the spatial class filtering for the whole 
study area after wheat was classified, Alphanumeric theme prints of 
classification results with and without employing spatial filtering were 
compared using aerial photographic mosaic as "ground infoiunation". Tlic 
combination of (2,2) was then selected as tho best t^hrcshold values in spatial 
filtering and was applied to five test sitps (c 40 km each) with different 
wheat densities for a quantitative evaluation ‘of tho accuracies of wheat area 
estimates. T tests showed tliat filtering With thresljold values (2,2) 
significantly decreased errors of commission and omission, also, tlie accuracy 
in area estimate was improved from the over estimate of 4.51 to 2.7^ and tho 
root-monn-squnro error dccroascd from 12C.18 ha to 107.02 ha‘. Extrapolating the 
same procedure of automatic classification using spatial filtering for post- 
classification to tho wliolo study aron,tlio accuracy in area estimate vms improved 
from the ovoresuimatc of 10.91 to 9.71. This study concltidcs _ that whdn a single 
pass LANDSAT data is used for crop identification and area estimation the post- 
classification procedure using a spatial filter provides a more accurate area 
estimate by reducing classification errors. 

1. INTRODUCTION 

The ideal LANDSAT pass for wheat identification and area estimate is in 
the late September or in tlic beginning of October when this crop is at yellow- 
ripe stage in soutliern Brazil. However, those months coincide w?th the initial 
of the rainy season, thus, cloud-covcr is ‘a serious problem for image analysis. 
LANDSAT data prior to the yellowing stage may also be used for area estimation 
purpose; ncvcrtlicloss ,aloss accurate estimate is obtained, since a large 
commission error is expected due to the similar spectral responses among wheat, 
pasturclnnd and fallow fields . The obj ctive of this study was to verify whether , 
improved classification results and a more accurate estimate of crop area might . be 
obtained through post-classification using spatial class filtering when only 
LANDSAT digital data of a single pass were available for analysis. 

2. STUDY AREA AND DATA ACQUISITION 

Cruz Alta is one of the major municipals for wheat production in Rio Grande 
do Sul State, Brazil. Tho geographic location of this municipal is around 28°35'S 
and 53°45'W. An area in Cruz Alta (approximately 20 x 40 km*), which represents 
tho wheat plantation of the state, was selected for this study (Figure 1). In 
this region, depending on climatic conditions, wheat may be planted in April or 
May and be harvested in October or November. 


•Presented at the Sixteenth International Symposium on Remote Sensing of 
Environment, Buenos Aires, Argentina, Juno 2-9, 1982, 
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Aircrnft data acquis ition 

On Soptombei^ 2, 19V9» INPE’s (Institiito do Pesquisns Efjpnciais) aircraft 
Bandoiraiito \^as flown ovor tho study area uslnfl a HC-10 photoBraimnetric camera 
and color xnfrarod (CIR) aoriul pbotoiu’nphs of medium scale (1:20,000) with 
30$ sldolappiUR and overlappinfi wore taken. Aerial photographs wore later 
interpreted visually according to a pj edotenninod legend. The 
photointorpretation results were then used as "ground information" to assess 
the computer-aided classification performance of Image- 100 system. 

LANDSAT data acmiisition 

Tho ideal hANDSAT pass for wheat identification and area estimation is in 
tho late of September or in tho beginning of October wlion wheat is matured 
and presents o golden-yellowish colors different from tho stirrounding crops 
(predominantly pasturoland) which are still green (Chon ct al, 1981). Ilowovor, 
LANDSAT data on September 22 , 1979 was with 1001 cloud cover, thus LANDSAT 
CCT's on September <ith when most of wheat plantations wore at hoading/floworing 
Stage, were used for this study. 

* . 3. METHODS ■ * . 

For wheat classification, the unsupervised clustering algorithm (K-moans) 
was first employed to separate homogeneous spectral classes; these spectral , 
classes Were then transformod to informational classes and training areas for 
onclt informational class woro located on the imago monitor of Image-100 system 
using tlio electronic cursor. Tho spectral information of those training areas 
wore used to derive training statistics Voquivod by MAXVlil^ wliich is a 
supervised classifier based on tho Gaussian mnximum-likoliliood decision * rule 
(Velasco ct al , 1978). This hybrid procedure of using clustering to assist the 
selection of training areas which were later used in a supervised 
classification was called M-2 procedure (Limn ct al, 1982). For homogonoination 
of classification results, n post-classification procedure named UNITOT was 
used 


UNITOT is a tliroo-by-throo pixel spatial class filter implanted in INPE's 
Imngo-100 system (Dutrn, 1982). There are two thrcsliold values wliiUi sliould bo 
predetermined by the analyst. The first threshold value, T-j , is the number of 
times tho analyst wants the central pixel to be considered in calculation of 
class frequency. After calculation of the .froquoncy for all the classes in the 
3x3 pixel matrix the highest class frequency will compared to n second 
threshold value Tj which is arbitrarily assigned by the analyst. If the 
value is smaller than the highest class frequency then the class of the central 
pixel will bo substituted by the class wlticli has the highest class frequency. I'f 
tho T, value is larger than the iiighcst frequency then the cl.ass of the central 
pixel remains unchanged. The best combination of Ti and T, to be used in spatial 
filtering should improve classification performance by diminishing 
classification errors. To select the optimal threshold values for this study 
the following combinations of (Ti , T,) wore tested, they were (1,2), (2,1), 
.(2,2), (2,3), (3,2), (3,4) and (3,5). 

Analysis procedure 

In order to work at the scale of 1:100,000 on the imago monitor th^e study 
area was divided into two square subareas (A and B) of approx. 20x20 km each. 
Analysis procedure was carried out similarly for both areas, Once the subarea 
was delimited, wheat was classified using the M-2 procedure as mentioned 
above, afterwards, spatial filtering (UNITOT) using different combinations of 
threshold values were applied and alphanumeric print out (1:20,000) of the 
classified wheat, with and witliout using spatial filtering, were obtained. 

Each alphanumeric printout v/as overlaid on the aerial photographic mosaic; on 



omsiKAi. la 

OF POOR quality ' 


ft Xiftht table i nml a visual cojnpnrison was iimdo obsorvinn tlio coiumissipu ftnd 
omission errors presented on the prinl-out. After compftring oil the print-outs 
with ftoriftl photogrttphic mosaic the host combination of throsliold viiluc was 
soloctcd. 

f 

In order to evulunto quantitatively the effects of spatinl class 
flltcrinR using the best threshold values on the clnssificotion results, five 
test sites 40 kiif eaclO with different whent densities wore selected from tho ‘ 
study area. A lioint-by-point comparis'on of nlphanumoric prlnt*outof each tost 
site to its corresponding aerial photor.rapliic mosaic provided dnta for 
statistical analysos. Paired ■t^^tests were applied to tho porcontaRos of 
correct classification COC) , error of commission ChC) and tho cstiiiuitcd wheat 
fti’oas obtained by using, and without using tho spatial filtering. Correct 
classification, commission errer) omission error and relative difference «re 
defined ns below t 


correct classification CCC I) n 


n9 of wheat pixels which wore classified 
corr e ctl y 

area-transformed pixel n9 of wheat from 
aerial pliotographs. 


commission error 


n9 of *non-whont pixels which wore erroneously 

(cn %) u 

n9 of wheat pixels classified by Imago-100 
system 


omission error (OH $) » 1 - CCl 


rolatlvo lUrfcrcnco (lilU) . w ..1)x10(ft 

estimated wheat area from -ncriol 
' ’ photographs 

4. lUlSUbTS AND DISCUSSION 

t 

Piguro ?, shows the clossificntion results of subarea A using procedure* 

M-2. liven thougli M-2 was selected as the best procedure for wheat 
clossification Chima at al, 1982> tliero were some confusion among tlio classes 
of wlicat, pasturclnnd and sparse arboreal vegetation. This might bo explained 
by the fact the l.ANDSAT data were gathered at t)ic time when wheat was in 
honding/Clowcring stage, consequently perfect sepn/ntion could not bo obtained 
due to the similarity between spectral responses of wheat and pasturoland. It 
was. noted that in most of tl\e cases the sparse arboreal vegetation which was 
misclnssified ns wheat possessed n understovo of vigorous grass whirh might 
contribute to the spectral similarity to wheat. 

Quantitative comparisons of nlphamimoric printouts to aerial photographic 
mosaic showed that the applicntion of spatial class filtering with the 
combination of tliroshold values (2,2) gave better classification results. This 
impYovement can bo observed by comparing Piguro 2 and 3. Quantitative 
comparisons wc ,'c performed using data of Table 1 and 2 which were obtained 
after classification with and whithout the application of UNITOT, combination 

(2.2) , for five to^<t oveas, T tests showed that the application of UNITOT 

(2.2) increased significantly (o r 0,05) tho porcontago of correct 
classification and decreased significantly (a « o.oi) the error of commission. 
Comparing area estimates obtained by using of without using UNITOT to that of 
tho aerial photographs no statistical difference was found. However, with tlio 
applying of UNITOT after classification the root-mcan-squaro error of area 
estimates for tltoso five areas decreased from 1 26. Ifi lift to 107,02 a. Tho relative 
difference in area of these five test sites decreased from n>< overestimate of 
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4.5$, when no spatial class fi'ltoi'ing was used, to +2.7$ when UNITOT with 
threshold values C2,2) was applied in analysis proccduro, Applyiny UNITOT 
(2,2) to the whole study area an improvomout in area ostinmtion was obsorvod, 
i.c. UD of +10.9$ obtained by classification whitliout the application of 
spatial class filforing was docroased to +9.7$. nven though the improvoments 
in area ostiiimto wore not pronounced the better classification accuracy caused 
by a liighor correct classification and a smaller ori'or of commission when 
spatial class filtering was employed made tjiis post-classification worthwhile. • 

5. CONCLUSIONS 

When multitcmporal LANOSAT digital data arc not available for crop 
identification studies the classification errors caused by the similarities of 
spectral responses among classes on a single-pass LANUSAT data may bo 
diminished by the applicatiOA of spatial class filtering with optimal threshold 
values in postvclassificntion. Tho improvoments in crop classification results 
using spatial class filtering for post-classification are contributed by 0 , 

.higher percentage of correct classification, a smuAlor error of commission and 
a more accurate estimate in area. 



Fig. 1 Location of study area. 
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FiR. 2 . Kctult of wheat c laisl f ient Ion without using UNITOT (2.2) in suharca A. 
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I-IS. S - KcsoJt of wheat cliiMi i f icat ion usiriii UMIOI (2,2) in nuharra A. 
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Table I, Ilstlinatotl wheat area based on aircraft and hANDSAT 
data for five tcsf aitoa 



TUST SITI; 

WlUAT 

Di^Nsrry 

Ci) 

1 HSTIMATHD V.’IIBAT AUHA (ha> 

1 

ArsIUAL 

PHOTOGRAPHS 

1 HANDSAT CCTs 

WITHOUT 
SPATIAh 
1-IhTKRING 1 

SPATIAh 
PlhTRRING 
WITH UNITOT 
(2,2) 

Bi 

so. 69 

2331.77 

2290.26 

(-1.781* 

2286.98 

(-1.92) 

2 

» 

46.13 

1897.60 

1906. '26 
(+0.451 

1852.51 

(-2.38) 

3 

39.28 

161S.72 

i 

1999.77 

(+0.51) 



1600.79 

(-0.92) 

4 

‘ 22.87 

971.05 * 

1 1481.95 

I ' (+23.551 

1137.96 

(+17.18) 

« 

5 

31.15 

1 

1322.50 

8502.32 
■ (+12.05) 

1481,06 

(+11.99) 

TOTAI. 

! 

i 

8138.64 

8502.32 

(+4.401 

8359.30 

(+2.7) 

WISH 


- 

126.18* 

107,02 

* Porcentogo of relative difference (RD'») , 


Table II. Wheat classification accuracies resulted by using and without 
using UNITOT (2,21 in five tost sites 


CORRHCT CLASSIPICATION 
iCC%) 

COMMISSION ERROR 

im) 

WITHOUT 

SPATIAh 

PIITORING 

■■ 

WITHOUT SPATIAh 
FlhTERING • 

SPATIAh 

PlhTERING WITH 
UNITOT (2,21 

83.75 

86.05 

14.73 

12.26 

82.83 

83.93 

17.54 

14.67 

77.51 

78.93 

22.89 

20,33 

88.41 

89.08 

28.45 

23.98 

87.85 

88.30 

21.60 

21.15 
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